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Abstract 

This work presents a rigorous framework based on coarse-graining to analyze highly compressible 
turbulence. We show how the sole requirement that viscous efltects on the dynamics of large-scale 
momentum and kinetic energy be negligible — an inviscid criterion — naturally leads to a density 
weighted coarse-graining of the velocity field. Such a coarse-graining method is already known in the 
literature as Favre filtering; however its use has been primarily motivated by appealing modelling 
properties rather than underlying physical considerations. We also prove that kinetic energy injection 
can be localized to the largest scales by proper stirring, and discuss the special case of buoyancy- 
driven flows subject to an external spatially-uniform acceleration field. We conclude that a range of 
scales can exist over which the mean kinetic energy budget is dominated by inertial processes and 
is immune from contributions due to molecular viscosity and external stirring. 
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1 Introduction 



This paper is the first in a series which studies the physical nature of compressible turbulence. The 
aim here is to give a systematic, theoretical approach based on coarse-graining (or filtering) to analyze 
non- linear scale interactions in such flows. It builds upon previous work of [lOj and |S1 H]. There 
are several motivations for this work. First, there is no unique way to specifying a notion of scale, 
such as defining large-scale momentum and large-scale kinetic energy, in compressible turbulence. The 
traditional approach in this subject has employed density-weighted averaging, also known as Favre 
averaging, to decompose a flow into large-scale and turbulent components. Using Favre averaging, 
density, /9(x), and velocity, u(x), are combined to yield a large-scale momentum, (pu), and a large- 
scale kinetic energy, i(pu)^/(p), with (...) denoting a space-average y Jydx(...). However, such a 
decomposition has been borne out of convenience to modelers and practitioners rather than physical 
considerations. It seems that a priori there is no fundamental reason to favor Favre averaging over any 
other combination of p and u, from an infinite number of possibilities. 

Second, and more importantly, this paper provides the foundation for subsequent work in which we 
attempt to answer basic questions about the physics of compressible turbulence. While the classical 
ideas of Richardson, Kolmogorov, and Onsager form the cornerstone for our modern understanding of 
incompressible turbulence, there are no grounds for extending such a theory to compressible flows. The 
potent ideas of an inertial range and universality are often invoked without physical basis in compressible 
turbulence. Elementary questions on the possible existence of a scale-range which is immune from direct 
effects of viscosity and large scale forcing, on whether energy is transferred to small scales through a 
cascade process, and whether such a cascade is local in scale, all remain unanswered. Resolving these 
questions is necessary to warrant the concept of an inertial range and to justify the existence of universal 
statistics of turbulent fluctuations. Furthermore, Kolmogorov's 4/5-th law for the energy flux is an exact 
result which has no counterpart in compressible turbulence. An analogous result would be essential for 
attempting to predict the scaling of spectra and structure functions. 

This work is also very useful from the standpoint of numerical modelling. Compressible flows, espe- 
cially in astrophysical systems, often involve a huge range of scales which cannot be simulated directly. 
The coarse-graining approach provides a theoretical basis for constructing models of turbulence that 
may faithfully reflect the dynamics at unresolved scales. The formalism that we employ is the same 
as that used in large-eddy simulation (LES) modelling of turbulent flows. This work thus provides a 
theoretical complement to those modelling efforts. However, while the equations we analyze coincide (to 
a considerable extent) with those that are employed in LES of compressible turbulence, their use here 
will be for rather different purposes. In LES, plausible but uncontrolled closures are adopted for the 
subscale terms, whereas the aim here and in subsequent papers, following |3l [4], is to develop several 
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exact estimates and some general physical understanding of these terms. Another difference is that LES 
generally takes the scale parameter ^ to be a fixed length of the order of the "integral scale" L. Our 
interest here is rather to probe all scales in the flow, including limits of small £ <ti L. 

In this paper, we shall show how a decomposition based on Favre filtering (and averaging) comes 
out naturally from the physical requirement that viscous effects have a negligible role in the dynamics 
of large-scale momentum and large-scale kinetic energy. We call this the inviscid criterion. Using the 
coarse-graining approach, we will prove rigorously the existence of an intermediate range of scales over 
which viscous dissipation and external kinetic energy injection can be made to vanish. Our decomposition 
also leads to two terms responsible for transferring kinetic energy across scales and constitute the so- 
called "subgrid scale flux," which we shall discuss in detail. 

The outline of this paper is as follows. In §[5] we present preliminary definitions and discussion. In 
§[3] we show how viscous dynamics can be isolated to the smallest scales by a proper scale-decomposition 
and in §0] we discuss the inertial dynamics based on such a decomposition. In §[5] we prove that it is 
possible to localize kinetic energy injection to the largest scales in a system and discuss buoyancy-driven 
(Rayleigh- Taylor) flows as a special case. In §[S1 we examine contributions from compressibility effects 
to the flux of kinetic energy across scales. We conclude with §[7] and two appendices, [X] and [B] which 
contain detailed proofs of our results. 

2 Preliminaries 

2.1 Governing dynamics 

In this paper we study the dynamics at various scales through a direct analysis of the compressible Navier 
Stokes equations without the use of any closure approximation. The equations are those of continuity 
([T]), momentum and either internal energy ^ or total energy 



Here, u is velocity, p is density, e is internal energy per unit mass, E — |up/2 + e is total energy per 
unit mass, P is pressure, /i is dynamic viscosity, F is an external acceleration field stirring the fluid, 
q = — kVT is the heat flux with a conduction coefficient k and temperature T. For convenience, we 
have assumed a zero bulk viscosity even though all our analysis applies to the more general case. We 
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have assumed that ^ = vp in independent of x. In general, /i(x) ^ r(x)", hence our assumption should 
hold in isothermal flows. It is also reasonable in flows which are turbulent enough such that temperature 
conduction is enhanced, making /i a weak function of x. 

2.2 Coarse- graining 

We first present a general approach to analyzing scale interactions in a turbulent flow. Following [TU] 
and |4j, we use a simple filtering technique common in the LES literature to resolve turbulent fields 
simultaneously in scale and in space. 

For any field a(x), a "coarse-grained" or (low-pass) filtered field, which contains modes at scales > ^, 
is defined in d-dimensions as 

a,(x) = j d'^r G,(r)a(x-t-r), (5) 

where G'(r) is a convolution kernel. It can be any real-valued function which is sufficiently smooth, decays 
sufficiently rapidly for large r, and is normalized so that j d'^r G{r) — 1. It is assumed furthermore that 
G is centered, / d'^r rG{r) ~ 0, and with the main support in a ball of unit radius, / d'^r |rpG(r) = 0(I). 
Its dilation in d-dimensions, G({r) = £^'^G{r/£), will share these properties except that its main support 
will be in a ball of radius £. If G(r) is also non-negative, then ^ may be interpreted as a local space 
average. Note that G(r) can be chosen so that both it and its Fourier transform G(k) are positive and 
infinitely differentiable, with G(k) also compactly supported inside a ball of radius I about the origin in 
Fourier space and with G(r) decaying faster than any power r~P as r — > od. See for instance Appendix 
A in [5] for explicit examples. It can be shown that for any kernel G(r) with the above properties, a 
coarse-grained function /^(x) is infinitely differentiable 0. 

We can also define a complimentary high-pass filter which retains only modes at scales < ^ by 

a^(x) = a(x) -af(x). (6) 

In the rest of our paper, we shall take the liberty of dropping subscript £ whenever there is no risk of 
ambiguity. 

It has been remarked by [17 and []Garnieretal09 that filtering a strong discontinuity in a field, such 
as an external shock from an explosion, introduces "parasitic" contributions which can overwhelm the 
turbulent fluctuations at small scales. Since our primary purpose in this work is a fundamental physical 
understanding, rather than modelling, of non-linear scale interactions, and since such a strong shock 
would interact with the flow, we consider it only natural to include its contributions to sub-scales < £. 

The filtering operation ([SJ is linear and commutes with space (and time) derivatives. We can apply 
it to the continuity and momentum equations ([IJ-© to describe dynamics of large-scale fields. However, 

^Under the very weak requirement that cix|/(x)| < oo over the domain Q of the flow. 
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there is no unique way to filter these equations. For example, we may define a large-scale momentum 
field either as 'pgUi or as pu^. Similarly, a large-scale kinetic energy may be defined as ^p^ju^p or 

3 Identifying the viscous range 

A key idea of this paper is that the scale-decomposition of momentum and kinetic energy should satisfy 
the inviscid criterion^ i.e. it should guarantee that viscous contributions are negligible at large enough 
length-scales. This is necessary for the study of inertial range dynamics if such a scale-range exists in 
compressible turbulence. 

3.1 Scale decomposition 

Under the assumption that /x(x) — fi is independent of spatial position x, we can coarse-grain eq.(l2]) 
and commute the filter with space derivatives in the viscous diffusion terms. We obtain 

_ 2 

dtWi + djipu^Uj) = -diP + ^l^j{[{^JU^ + diUj) - -9™?!™%]} + pF^. (7) 

With such a decomposition, the functional form of viscous terms in ([T]) is similar to their counterpart in 
incompressible flows. If /(ix|u(x)p < cx), it can be shown rigorously that each of the viscous terms in 
eq. ([7]) is bounded by (const. at every point x. Therefore, the large-scale momentum, defined as 
pu^, does not diffuse under the action of molecular viscosity when /i/^^ <C 1. The proof is standard in 
real analysis and for applications in turbulence theory, see [5] and [1] . We repeat the proof in Appendix 
A for completeness. The idea behind it is simple; a term pV^u may be expressed in terms of a big-0 
bound, 0(/i (5u(i')/^^), which becomes negligible as p —> 0. Here, an increment is 5u(.£) = u(x + £) — u(x). 

The reader might question the value of a careful proof when one can arrive at the same conclusion 
by a simple dimensional argument. To illustrate the potential pitfalls of dimensional reasoning here, 
consider the quantity pVu : Vu^ . It may be argued that this should also vanish as — ?> for some fixed 
£ > 0. However, it is well known in turbulence literature that it does not (see for example [20l fT6] ). 
The problem lies in that /iVu : Vu£ cannot be expressed in terms of quantities at scale such as 
0{p5u{tY / 1"^), as one might innocently expect. When p — > gradients can become unbounded and the 
term Vu : Vu£ diverges. Even though Vu : Vu^ has small wavenumber modes < K ^ ^ it can be 
dominated by contributions from u(q) with wavenumbers Q ^ K due to the convolution Vu : Vu in 
Fourier space. Put more explicitly, while the product Vu(q) : Vu(k — q^Q^v^^^i^-q.)-^ j^^g g, Fourier 
mode at wavenumber |k| < K, it is proportional to ~ Q^. This example has a direct bearing on our 
definition of large-scale momentum. If we were to define large-scale momentum as p^u^ rather than 
pu^ as we did above, a viscous term in the balance equation would have the form /ip^p^^V^u^. Here, 
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again, the filtering operation would not commute with the laplacian and, due to possibly dominant 
contributions from high wavenumber modes ^ we would not be able to guarantee a priori that 
viscous terms are negligible at large 

In the analysis above we needed to assume that /i is a constant. The results should still hold 
if /i(x) only varies over large length-scales ^ £ or if its fluctuations are small relative to the mean 
|^(x) — ^ (/i). A priori tests of compressible turbulence by []Vremanetal95, [jVremanThesis, and 
[jMartinetalOO seem to suggest that, indeed, the additional non-linearity in viscous terms introduced by 
Sutherland's law is small. 

The scale decompostion employed in the large-scale momentum balance ([7]) is equivalent to traditional 
Favre filtering (see for example [9J), where a Favre filtered function is weighted by the density: 

/.(X) ^ (8) 

The operator (~) is linear but does not commute with derivatives. The large-scale momentum balance 
([7]) can be rewritten using definition ([8|) as 

dipui -f djifiu., Uj) = -dj (p T{ui,Uj)) - diP 

2 ~ 
+ ^l^j{[{^jUi + diUj) - -dmUmSij]} -I- 'pFi. (9) 

This is the same as the "bare" momentum equation ^ itself but with an additional contribution from 
turbulent stress, 

pT{Ui, Uj) = p{lL[Uj - Ui Ui), (10) 

which accounts for the effect of eliminated scales < I and vanishes identically in the absence of fluctua- 
tions at those small scales. One can also obtain a continuity equation for large-scale density: 

9tP + a,(pM0 = 0. (11) 

A main advantage of the filtering approach to analyzing turbulent fiows is an ability to resolve the 
relevant physical processes both in scale and in space as is apparent from the balance eqs. (|9|). (|lip . 
They describe the evolution of large-scale momentum and large-scale density at every x in the Sow and 
at variable resolution I. Using p£ and pu^, and eqs. ([9]) . (fTT|) . it is also straightforward to derive a simple 
dynamical equation for the kinetic energy density at scales > ^, for arbitrary L This yields 



dt-p}^ + V-J, = -n, ~ A, + PiV-n, -De + e^^ (12) 



where Je{x) is space transport of large-scale kinetic energy, + A^ (x), which we examine closely in 

§[31 is usually called the subgrid scale (SGS) kinetic energy flux to scales < ^, — P^V-u^ is large-scale 
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pressure dilatation, Diipc) is viscous dissipation acting on scales > £, and e™"'(x) is the energy injected 
due to external stirring. These terms are defined as 



n,(x) 

A,(x) 

J,(x) 
el'^ (x) = Ui pFi 



p djU, T(Ui, Uj) 



= A* 



p—^Uj + Fuj + UtpT[Ui, Uj) - puidjUi - -UjdiUi 



(13) 
(14) 

(15) 

(16) 
(17) 



where we employed in ()14p the notation 

T,{f,g)=JM',~J,g, (18) 

for 2"'^-order generalized central moments of any fields /(x), ^(x) (see |10]). 

Just as we have shown that viscous diffusion of large-scale momentum pu = pu is negligible, we can 
also rigorously prove under very weak conditions that viscous dissipation De(x.) of large-scale kinetic 
energy ip|up vanishes at every point x when p/£'^ ^ 1. The proof, given in Appendix A, uses standard 
tools from real analysis. 

3.2 The inviscid criterion 

In one of his original articles, [7] motivated the usage of density-weighted averaging by the fact that 
average mass of a fluid in a volume V advected by the large-scale velocity (pu)/(p) is conserved. We 
shall now briefly repeat Favre's argument. The change of average mass in a volume V advected with 
some large-scale velocity u* is 

/ d^^ dt{p) + V- i{p)u*) = I d^x V- ((p)u* - (pu)) , (19) 
Jv Jv 

where now (...) denotes ensemble averaging. The equality follows from using the ensemble averaged 

continuity eq. ([T]). Favre averaged velocity u is defined as the choice of u* which makes average fiux of 

mass due to fiuctuations (or turbulence) vanish, (p)u — (pu) 0. The same argument carries over to 

spatially filtered dynamics, where now u^ is defined as the choice of large-scale velocity which does not 

lead to subgrid terms in the filtered continuity eq. dTTI) . 
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There are two comments wc would like to make concerning Favre's argument. First, the special 
property large-scale velocity u enjoys, i.e. suppressing turbulent diffusion of mass in eq. (|lip. does not 
logically imply by itself that a density-weighted decomposition is a necessary choice in the mass balance 
([TT|l . It is certainly not unphysical for a turbulent flow to diffuse mass and, in this respect, the choice of a 
large-scale velocity field would depend on the particular purpose of an investigation. Second, a criterion 
requiring that turbulent mass diffusion be zero has no logical implication on the scale decomposition of 
momentum and kinetic energy. Unlike in eq. (Ilip . the Favre decomposition results in turbulent diffusion 
and dissipation of large-scale momentum and kinetic energy as seen from eqs. (P|). (|12l) . 

Yet, density- weighted filtering is used extensively in LES of compressible turbulence due to its mod- 
elling appeal. One of the perceived advantages is the absence of subgrid terms to be modelled in the 
coarse-grained continuity eq. PT|) . Another reason is that Favre filtered equations (p]). ([T^ are struc- 
turally similar to their classically filtered counterparts in incompressible flows, which allows for carrying 
over models from the incompressible LES literature. Furthermore, none of the subgrid scale terms 
is a function of pressure which practitioners try to avoid modelling. When using the ideal gas law, 
P = (const. )pT, as the equation of state, there is also an added advantage that filtered pressure, P, can 
be expressed as a function of resolved quantities, p and T, without additional subgrid terms. 

What we have shown above is that Favre decomposition of momentum and kinetic energy satisfies 
the inviscid criterion. It guarantees that viscous contributions are negligible at large enough length- 
scales. Such a decomposition of momentum and kinetic energy is borne out of a physical requirement, 
irrespective of practical modelling considerations. We remark, however, that it may not be the unique 
decomposition satisfying the inviscid criterion. In other words, we did not prove that it is necessary. We 
only showed that the Favre decomposition is sufficient to satisfy the inviscid criterion. 

As we mentioned in the introduction, while our equations (l9|). (fTT|) . (|12p coincide to a considerable 
extent with those that are employed in LES of compressible turbulence, their use here and in the 
ensuing papers will be for rather different purposes. Whereas the primary goal in LES is to model the 
subgrid terms, the aim here is to develop a physical understanding of these terms and estimate their 
contributions at different scales, including limits of small £ <ti L, through exact mathematical analysis. 

4 Inert ial dynamics 

Now that we have isolated viscous effects to the smallest scales where is defined as the scale at 
which viscous effects become significant in kinetic energy balance (|12p. we can study the dynamics at 
scales £:s> 
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4.1 Deformation work 

The first term in kinetic energy SGS flux, Hi in eq. (jl2p . is the contribution from deformation work 
done by large-scale strain djUi against the subgrid stress 'pr^Ui, Uj) (see for example |22|). This is similar 
to its incompressible counterpart except that the strain is not traceless here. It acts as a sink in the 
large-scale kinetic energy budget ([T2l) and represents that part of the kinetic energy transferred from 
scales larger than £ to smaller scales at point x in the flow. 

Furthermore, is Galilean invariant due to the subtracted large-scale terms in definition (jlOp of 

the turbulent stress. Other definitions of the SGS fiux are possible such as Uidj{pT{ui,Uj)) which differs 
from our definition ()13|) by a total gradient dj{'pUiT{ui,Uj)). However, this definition is not pointwise 
Galilean invariant, so the amount of "energy cascade" at any point x in the fiuid according to such a 
definition would depend on the observer's velocity! [TT], [TB], and [TU] all emphasized the importance of 
Galilean invariance. More recently, |6] and [2] showed that Galilean invariance was necessary for scale- 
locality of the cascade. There are non-Galilean-invariant terms in our budget (|12p but, as is physically 
natural, they are all associated with space transport J of kinetic energy. 

Another physical requirement on the flux Hi (x) is that it should vanish in the absence of fluctuations 
at scales smaller than £ (or a moderate fraction thereof); for example, when £ is equal to 2TT/Kmax, 
where Kmax is the maximum wavenumber in a direct numerical simulation (DNS). This is satisfied by 
our definition of (x) identically at every point x in the flow. Other definitions of an energy fiux are 
possible, such as the "unsubtracted flux" of an incompressible flow critiqued in [6] and [2], 

n^"" (x) = u, JTC, , J\fC, = dj {pu^uj ) 

which is often employed in literature that considers the sharp-spectral filter. Here, AfCi denotes the 
nonlinearity in the momentum equation. 

Using such a filter, the "unsubtracted flux" across wavenumber K is computed as 

n^'-(x) = Yl MpV"' E A^«(k)e^'^^ 

\p\<K \W\<K 

Taking K — Kmax, we have H^^^(x) = Uidj{puiUj) which is in general nonzero. It is only after 
averaging over all space (or in the absence of flow beyond the domain boundary) that one gets 

(nFlJ = (5,(p^^,))= E ^:(k)A^.(k)-o. 

Similar considerations apply for compressible flows, where an unsubtracted flux may be deflned as 

Hr^(x) = u,JTC,~ AT, UC, - d,{pu,u,), M = d,{pu,), 

with ^>fCi and A/" denoting the nonlinearities in momentum and density equations, respectively. This 
flux does not vanish in general when £ = i^'^a^,, except after space-averaging. 
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4.2 Baropycnal work 

The other part of kinetic energy flux, Kg in eq. (jl2l) , is intrinsically due to compressibility effects and 
vanishes in the absence of density variations. It represents work done by a large-scale pressure- gradient 
forct^ 'p~^'^ P against suhscale mass flu^T{p^ u). We shall refer to A£(x) as baropycnal work due to its 
inherent dependence on pressure and density variations. It is not entirely due to baroclinic effects for 
it can be non-zero even when small-scale density variations and VP^ are aligned as shown in Figure [1] 
Similar to II^, it also acts as a sink in the large-scale kinetic energy budget (IT^ . is pointwise Galilean 
invariant, and vanishes identically at every x in the absence of fluctuations at scales < i. Baropycnal 
work is known to play a major role in turbulent combustion (see for example |21| and jl2|). It has also 
been recently observed by |13) and |14) to play a major role in turbulence production in buoyancy-driven 
flows with significant density differences. In such flows, termed "variable-density flows" by |13j . there 
are two or more incompressible miscible fluid species, such as water and brine, which have significantly 
different densities. 

The physical mechanism behind this part of the SGS flux, illustrated in Figure [TJ is simple. In a 
frame co-moving with a ball of radius t in the flow, a pressure-gradient force 'p^^'^Pi at scales > £ 
acts on the ball within which the fluid has non-uniform density (density variations at scales < t). Per 
unit time, if p2 > pi in Figure [T] and |ui| — |u2| = u, then parcel 2 gains 'p~^'V Pep2U in kinetic energy 
from the large-scales and parcel 1 loses 'p^^'V Pipiu to the large-scales. On aggregate, scales < i in this 
ball would gain p^^VP^ (p2 — Pi)u in kinetic energy from the large-scales. This effect would vanish in 
the absence of density variations. Note that if |ui| does not equal |u2|, then the whole ball will have 
an average large-scale velocity = Ui + U2 which does not play a role in such a process of inter-scale 
energy transfer — hence the relevance of the premise of a co-moving frame. 

Similar to deformation work, 11^ (x), baropycnal work, A£(x), is not sign-deflnite in space. We expect 
that at points x in the flow where large-scale pressure gradient VP£(x) opposes the density gradient 
due to small-scale fluctuations (p2 > Pi in Figure [T]), then a motion similar to that of a Rayleigh- 
Taylor instability would ensue such that velocities Ui and U2 will be as illustrated in Figure [TJ We, 
therefore, expect that at such points in the flow A£(x) will be positive and transfer energy to small 
scales. Whether, on average, the subscale mass flux t(/3, u) would correlate positively or negatively with 
large-scale acceleration fleld VP can be determined empirically. 

Despite the recognition of A^'s importance in turbulent combustion and variable-density flows, to the 
best of our knowledge, the presence of baropycnal work has not been appreciated in literature on LES of 
compressible turbulence. In the LES community, this term has always been lumped with P^ V-u^ in the 



^The term "pressure-gradient force" is used in the meteorology literature. It is not a force but an acceleration. 
^Here, "flux" denotes a flux in space, not in scale. 
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Figure 1: Heuristic explanation of the physics behind baropycnal work, A^. In a ball of radius there 
are two small-scale fluid parcels of equal size with densities pi and p2, and antiparallel velocities Ui and 
U2, respectively. In a frame co- moving with the ball, we must have |ui| = |u2|. 

form PfV'U^ (plus an additional space-transport term) and treated as a large-scale pressure dilatation 
which does not require modelling. We do not make any claim about the importance of baropycnal work 
relative to deformation work, which we believe is an issue that deserves careful empirical investigation. 
The reason we keep baropycnal work separate from pressure dilatation is due to a fundamental distinction 
between the two. Notice that both deformation work, Hi, and baropycnal work, A^, involve large-scale 
fields acting against small-scale fluctuations. This makes them capable of transferring energy across 
scales. On the other hand, pressure dilatation, Pi'V-ug, involves only large-scales and cannot transfer 
energy directly across scales. 

4.3 Pressure dilatation 

The pressure dilatation — PV-u in eq. (|12p represents conversion of large-scale kinetic energy to internal 
energy through compression. In the incompressible limit, this vanishes at every point x. Unlike 11^ and 
Ai, pressure dilatation only contains scales > £ (at least for filter kernels G(k) compact in Fourier space). 
Therefore, it is not involved in the transfer of energy across scale i and it does not vanish in the absence 
of subscale fluctuations. 

5 Kinetic energy injection 

Under statistically steady-state conditions and after space averaging, the large-scale kinetic energy bud- 
get (fT2|) at scales £^ £^ becomes 

(n,) + (A,) - (P.V-u,) = {eT'). (20) 
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We assumed that none of the kinetic energy is transported beyond the domain boundary. We also 
dropped viscous dissipation, (-D^), which we have proved to be neghgible. If, furthermore, the injection 
of kinetic energy can be made localized to the very largest scales i ^ £, then one should be able to 
study non-linear inertial dynamics at intermediate scales L £ ^ over which the left-hand side of 
eq. ((20|) becomes a constant, independent of £. 

The injection can be made localized to small wavenumbers if the external acceleration field F(x) 
has only small Fourier wavenumbers. To show this, we will denote a field in a periodic domain [0, 27r)'^, 
coarse-grained with the sharp-spectral filter to retain only Fourier modes |k| < K, by 

a<^(x) = ^ dk a(k)e*-^. (21) 

\k\<K 

This is similar to a£(x) with £ being of the same order as K~^. We will also employ in our argument the 
notation [Ki,K2] for a band of Fourier modes Ki < |k| < K2 and denote the corresponding band-pass 
filtered field by a\^^'^^^. Now consider an imposed acceleration F^^" with modes restricted to small 
wavenumbers < Kq <^ K. The net energy injection from (1171) is, thus, 

We have the following exact identity: 

{pF<^°)^^ (x) = (p<-^+^" p<Ko'^<^ ^ p<K~Kop<Ko _,_ ^p[K-KoM+Ko]p<KoY^ 

= p<'^F<^° - p^^^-^^°^^\ F<^° + (p[K-Ko,K+Ko\p<KoY^ . (23) 



The first equality follows from the fact that a Fourier mode /5(p) with |p| > K + Kq along with F(q) 
with |q| < Kq cannot contribute to the product pF with Fourier modes k = p + q restricted to |k| < K. 
The second equality is due to {^p<K-Kop<Ko'^<K ^ p<K-Kop<Ko^ Notice that in the second and third 
terms of the last expression in (j23l) , density has only modes in a thin band of width at most 2Kf) around 
K. As Kq/ K -H- 0, a factor p[K-KQ,K+Ka\ gj^ould become negligible due to a decaying density spectrum 
and, therefore, the second and third terms in ((23|) vanish. Using this observation in p2|) implies that 
the net energy injection with Kq <C K reduces to 

(^"') « {^^^^P^'^F^''") - ( (P^^)^""" F'"'")^ (24) 

where the last equality follows from orthogonality of Fourier modes. In Appendix A, we present a 
rigorous proof of the above statement under precise (and weak) conditions. Our result ([M)) implies 
the non-trivial possibility to make injection localized to wavenumbers < Kq by employing an external 
acceleration field band-limited to wavenumbers < Kq. 
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In arriving at eq. (j24p . we had to stir with an external acceleration field such that the force in 
is weighted by instantaneous density p(x,t). Had we stirred the momentum equation with an 
external force ^<^° instead of pF'^^° , the injection would have had the form (e"'-') = ( ^^^<k — •^i'^") ~ 
(^u^^o.^<Ka'^^ for Kq ^ K. We shall not present a detailed proof of the latter statement but the idea 
d details behind it are similar to those leading to (|24l) . Finally, we note that the acceleration field F^-'^" 
in momentum eq. ([2|) is of a general form and can have both solenoidal and irrotational components. 

Under such conditions, one obtains that the sum of kinetic energy flux and pressure dilatation, 

(n^) + (A^) - (PiV-ue) = (const.) (25) 

in (P0|) . is constant, independent of wavenumber K = £~^. This in itself is not an analogue to Kol- 
mogorov's 4/5-th law since it contains the pressure dilatation term which does not involve energy trans- 
fer across scales. However, we shall show in a subsequent paper that it is possible to derive a result 
analogous to the 4/5-th law under a reasonable assumption on the co-spectrum of pressure dilatation. 

5.1 Special case: Rayleigh- Taylor flow 

We shall apply our result on localized injection to the case of a Rayleigh- Taylor flow driven by grav- 
itational forces. We ask the following question: At what scales does mean conversion of gravitational 
potential energy to kinetic energy take place? 

We shall show that, in the presence of a spatially uniform gravitational field g(x) = g, such conversion 
into kinetic energy only takes place at the largest scale — that of the domain size. Consider the large- 
scale kinetic energy equation (I12p where F(x) is replaced with g in e™-*. We have from ^17} that net 
input of kinetic energy due to gravity is 

(e"'^) = {urPeie) = iWrg.) = g-(pu) (26) 

where the last equality follows from J dxfg{x) — J drGe{r) J dx/(x -I- r) = (/). Result (|26|) shows that 
mean injection is independent of scale £ and only takes place at the scale of the domain, Ldom- 

Put in more detail, the mean kinetic energy (p^ju^p) /2 at scales > £ increases at a rate (u^-p^gf) due 
to gravitational forces. Consider a sequence of scales li > £2 > ■ ■ ■ > £n > ■ ■ ■ and the average amount of 
potential energy being converted into mean kinetic energy, (e™"') = (u£„'P£^g£„), at successively larger 
sets of scales [£n, Ldom)- In general, as n — > 00, (e"'"') approaches the total amount of potential energy 
converted into kinetic form. The fact that (e™"*) is independent of ^„ over the entire scale-range [0, Ldom) 
implies that all conversion takes place at the domain scale. 

Our argument demonstrates the power of the filtering approaclj^ in analyzing non-linear scale inter- 

^ In this particular case, with a constant g containing only a fc = mode, the same conclusion would have also been 
possible by examining the injection into mean subscale kinetic energy, (p|up^ — p^lu^ P)/2, which may be easily determined 
to equal zero at any £. 
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actions. We were able to arrive at our answer precisely because the filtering technique allows for probing 
a continuous range of scales, in contrast to the traditional averaging approach. 

The conclusion is probably non-intuitive at face value because, in a Rayleigh- Taylor flow, "fingers" 
of heavy fluid at any scale penetrate the lighter fluid as they descend, converting potential to kinetic 
energy. It seems to contradict our result that conversion only takes place at the domain scale. Key 
to understanding such an ostensible paradox is the definition of kinetic energy based on a Favre scale- 
decomposition. Mean kinetic energy at scale i may be rewritten as (|pu^p/p^)/2, which emphasizes the 
central role of momentum in defining scale. Indeed, result (j26p demonstrates that it is only the A: = 
mode of momentum which participates in converting potential energy into kinetic energy at scale with 
mode k — 0, (pu)^/ {p)/2. While momentum at scale £ can have contributions from density and velocity 
at all scales — for example, mean momentum (pu) = ^j^. /5(k)u(— k) — the scale £ of kinetic energy 
depends on that of momentum pu^ and not on, for instance, pgUi. 

6 Compressibility effects in the SGS flux 

Our scale-decomposition allowed us to identify two SGS sinks for the large-scale kinetic energy budget, 
namely deformation work, 11^, and baropycnal work, Ai. 

While Hi also represents a cascade mechanism in incompressible turbulence, is solely due to 
density variations and thus is intrinsic to compressible fiows. However, even deformation work has con- 
tributions from compressibility effects. We shall show in a subsequent paper that 11^ — 'pdjUiT{ui, Uj) + 
. . . 8 terms . . . , using exact identities. If we consider that part of 11^ equal to 'pdjUif{ui, uj), we see that 
compressive modes with k-u(k) ^ can play an important role in deformation work. This is because 
the large-scale strain Vu is not traceless. The effect of compressive modes is best illustrated by the 
one-dimensional Burger's flow, 

dtu + dx{ — ) = vdxxU. (27) 
A large-scale kinetic energy budget analogous to can be derived (see for example [S]): 

d}^ + dx^^+uT{u,u)-vdx^-'^=-Tl''''''' -u\dxu\\ (28) 
in which viscous dissipation on the right-hand-side is negligible and 

Y]bnrg ^ __Q^- -(^^ ^) (29) 

is the only sink. There are no shearing motions in Sow (l?f)) . The only way energy cascades to small 
scales is through n^"''^ which is solely due to compressive modes. As is well-known, this cascade is 
manifested in the formation of shocks. 
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In general, a simple measure that quantifies kinematic role of compressibility on deformation work 
(USD is 

^comp = n, + p d,lF, T{ut, uj), (30) 

where the velocity u = u" + u"^ is decomposed into solenoidal and irrotational components, and u'^, 
respectively. 

7 Summary 

In this paper we have shown that viscous diffusion and dissipation can be isolated to the smallest scales 
in a compressible flow by using a proper scale-decomposition. Guided by this physical requirement, 
which we call the "inviscid criterion" , we found that a Favre decomposition of the momentum and 
kinetic energy into large-scale and small-scale components is sufficient to guarantee a negligible role of 
molecular viscosity in the large-scale dynamics of high Reynolds number flows. 

We were also able to establish through an exact analysis that mean kinetic energy injection can also 
be made localized to the largest scales in a flow by proper stirring. We discussed the special case of 
buoyancy-driven flows in which stirring is due to an external spatially-uniform acceleration held, and 
showed that mean injection of kinetic energy occurs only at the very largest scale in the system. 

Localizing viscous dissipation to the smallest scales and energy injection to the largest scales is 
necessary to allow for studying inertial dynamics at intermediate scales. Under steady-state conditions, 
satisfying these two ingredients implies that the sum of mean SGS kinetic energy flux and pressure 
dilatation, (II^ + Ki) — {VPi-ui), is constant, independent of scale £ over the intermediate range L ^ 
i ^ ifj,- This is simply a consequence of energy conservation; whatever energy is input into the system 
has to either reach dissipation scales by way of the SGS flux or get converted into internal energy through 
pressure dilatation. 

The results of this paper lay the groundwork for ensuing work to tackle basic issues concerning the 
physics of compressible turbulence. Questions on whether energy reaches dissipation scales through a 
local cascade process, on whether an analogue to Kolmogorov's 4/5-th law exists, and on how spectra 
and structure functions should scale in compressible turbulence, will form the subject of subsequent 
works which build on the conclusions of this paper. 
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A Viscous effects 

The following proposition shows that viscous terms in the balance ^ of large-scale momentum p^Uf are 
negligible when ^/£'^ <^ 1 . The proof is identical to that given in [5] and [1] . 

Proposition 1. If velocity solutions of the compressible Navier-Stokes equation {2p-([3^ over a domain 
fl have finite 2nd-order moments: Jj^dxju^p < oo, then viscous terms in the large-scale momentum eq. 
^ vanish pointwise as /i — > 0. 



Proof of Proposition 1: The viscous terms in eq. ([7]) can be bounded using Cauchy-Schwartz inequahty, 

A*V2lIF,(x) = ^ y"dr(V2G),(r)u^(x + r) 



<fl|(V^G).|y|u-(r)||^, 

with II = {\...\Py/P the L p-norm. Since G(r) € , its derivatives are (uniformly) bounded and 

for any fixed £ > 0, ||(V^G)£||2 < (const.). Hence, the viscous term ^V^u^^ -> at every point in space 
in the limit of vanishing viscosity. The same argument applies to /iW-u^. □ 



The following proposition proves that viscous dissipation of large-scale kinetic energy ip^ju^p is 
negligible when /i/^^ ^ 1. 

Proposition 2. If solutions (/3^,u^) of the compressible Navier-Stokes equation over a bounded 

domain Q. have finite 3rd-order moments: J^dx.\p^^\^ < oo and Jj^(ix|u''p < oo, then viscous terms in 
the large-scale kinetic energy budget I112\) vanish pointwise as fi ^ 0. 

Proof of Proposition 2: In the following, we shall drop superscript /i from (p'',u''). Using the exact 
identity u = u + r (p, u) /p, we have 

|9jMi(x)| < \d.fui \ + \^dff(^p,Ui)\ + \Xfr{p,Ui)djp\. (31) 

The first term \dfui\ is bounded by ^~^(const.)|| VG^ ||2||u||3 through an argument identical to that 
in Proposition [T] and using the fact that ||u||2 < (const. )||u||3 over a bounded domain. 
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The second term in (PT|) can be rewritten as 



J dr{djG)i{r)p{x + r)u,{x + r) 



dr(a,G),(r)p(x + r) J dr'G,(r')u.(x + r') - j drG,(r)p(x + r) J dr'(9,G),(r'K(x + r') 
using integration by parts. Employing the 3-3-3 Holder's inequality, this expression is bounded by 

rii||(VG),||3(l + 2||l|a|G,||3)||p||3||u||3 
The third term in ([5T|) can be rewritten as 



-r{p,u,)d,p=-p-'r' I dr(a,G),(r)p(x + r) 



J drGi>{r)p{y: + r)u,{yi + r) - J drGi{r)p{yi + r) J dr'G^ (r')ui(x + r') 
Using Holder's inequality, this is bounded by 



^"'^l|(VG).||3l|l|l3l|G43(l + l|l|l3l|G43)|IHl3ll"ll3 



Finally, the viscous term is bounded by 



fildjUidjuA < 



Hi 

12 \ 



(const.) + (const.) ^_ ^ + (const.)— 
P P 



The coarse-grained density Pf(x), for positive filter kernels G(r) > 0, is proportional to the mass in a 
ball of radius £ centered around x. At any point in the flow, this is bounded away from zero. Hence, 
for any fixed ^ > 0, the viscous term fidjUidjUi — > at every point in space in the limit of vanishing 
viscosity. The same argument applies to pdjUjdiUi. □ 

B Kinetic energy injection 

The following proposition proves that mean kinetic energy injection is localized to the largest scales 
> L. In the text, we presented an informal argument using the sharp-spectral filter. However, this 
filter's kernel is not positive in physical space, which prevents us from inferring that (p^^)^^(x) is 
bounded at every point x @. For this reason, we will use in our proof below the Fejer summation of 
Fourier series. This is defined as 

aW(x)^lf]a<^(x), (32) 



®We are not claiming that coarse-graining with a sharp-spectral filter yields non-localized injection, but only stating 
that we cannot prove otherwise. 
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where a<^(x) is coarse-grained using the sharp-spectral filter as defined in (l?T|) . The Fejer kernel r/<-(r) 
is positive in physical space and compactly supported in a ball of radius K in Fourier space. 

To show that mean kinetic energy injection is localized to the largest scales > -f^cT^, we need to prove 
that 

(4"^) - ^ ( (ipu^r^Y"' F<-')^ (33) 

as Kq/K 0. 

To this end, we will need an identity similar to (1^^ in the text for the Fejer filtered quantity. For 
K > 2Kq > 2 (the special case of K > 2Kq = is straightforward and was discussed in ij lS.ip . we have: 

-, K 2Ko-l 
Q=2Ko Q=l 

Using identity (|23p. the first term on the right hand side can be rewritten as 



K K 



K ^ ^' ' K 

Q=2Ko Q=2Ko 



2Ko-l , K 



= p(^)i^<^o _ J_ ^ p<Qp<Ko ^ i_ ^ „p[Q-Ko,Q]p<Ko _^ f^plQ-Ko,Q+Ko]p<Ko-^<Q^ 
^ Q=l ^ Q=2Ko 

We finally have the identity 



2Ko-\ 

\<Q _ <Qp<Ko 



Q=l 

L ^ _p[Q-Ka,Q]p<Ko (^p[Q~Ka,Q+Ko] p<Ka-^<Q I 
^ Q=2Ko J 



Proposition 3. If solutions (p, u) of the compressible Navier-Stokes equation (1^-^ over a periodic 
domain T ~ [0, 27r)^ satisfy pl^g = /^(ix|/3p < oo and (pu)^^^ = J,j,(ix|/9Up < oo, then 



hm l(e"'^- ((«w,:)('^)r'^V<^«) 



< Prras{pVL)rms 1™ | (cOnst.) W ^ -|- (cOnSt.) ^ i- = 0. (35) 



Proof of Proposition 3: 

We first note that prms < oo implies that the density spectrum, EP(k) = X)fc-o 5<|k|<fc-i-o 5 IpI'^)!^; 
is bounded by a power-law: 

EP{k) < (const. )p^„^fc"", n > 1. 
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Otherwise, pmis increases without bound as a function of Reynolds number (see for example [5]). 
Using identity we have 



{^k') - ( {{pud 



?<Ko 



/ {pUj)^^ I 1 ^ (^pp<Ko-^<Q_p<Qp<Ko^]_ ^ _^[Q-Ko,Q]^<-R'o_^(^[Q-ifo,Q+i^o]^<-R'o) 

Applying the Holder inequality to the right hand side, we get 



2Ko-l 



< 



(36) 



Since \\p^'^ Ko,Q+Ko] ||2 — 2 dkEP{k), we obtain from our bound on the density spectrum that 



< (const. )/9rms\/-f^0 



^ q' ^ q' 
1 



1/2 



This yields that 



}_ < (const.)p™./^-|X:Q 



i/2 



Q=2Ko 



< (const.) PrmsVKo 



dQ Q- 



V2 



KJ, 

(const.) Prms'\J~~^ ( ^(n-l)/2 ~ 1/2 ^ 

(const.)p™sy^(i^^ 



n 7^ 2 
n = 2 



(37) 



Since p^^^ is proportional to the mass in a ball of radius ^ J^^^^, we have || (p*-^-*) Hoc < oo. 
We can also have IjF^-'^'' ||oo < oo by construction, and ||(pu)*^^'||2 < (pu)r,ns < oo follows from our 
assumptions. We finally have from that 



K 



(pu. 



<K„ 



?<Ko 



)l 



< Prms (pu)rms< (cOnSt.)A / — + (cOUSt.) — 

which vanishes as Kq/K 0. □ 
REMARK: In our proof, we only used that the density spectrum is bounded by (const. )p^„jg/c^". 
This is very weak and does not even require a continuous spectrum. If the spectrum decays faster than 
such a bound, then (e*^^) ( ((pw,)^^^) ° F<'^°) at a rate faster than what is stated in estimate 
(El). 
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